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i!S ■ ABSTRACT 

o : 

^ , We present radial velocity and photometry for four early-type, massive 

^ \ double-hned spectroscopic binaries in the R136 cluster. Three of these systems 

j> \ are eclipsing, allowing orbital inclinations to be determined. One of these sys- 

\ tems, R136-38 (03 V + 06 V), has one of the highest masses ever measured, 

\ for the primary. Comparison of our masses with those derived from stan- 

■ ■ ■ ' dard evolutionary tracks shows excellent agreement. We also identify five other 
light variables in the R136 cluster which are worthy of follow-up study. 
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1. Introduction 

Empirical checks on the mass-luminosity relation and on how well the evolutionary 
models match reality are sorely lacking for high mass (> 20A4q) stars, even near the zero- 
age main-sequence (ZAMS). It has long been recognized that evolutionary tracks (past the 
ZAMS) are affected by mass-loss (see, for example, Brunish & Truran 1982), and it is now 
well recognized that in particular the models are highly sensitive to how mixing of material to 
and from the core is treated (Maeder & Conti 1994), with the most recent work emphasizing 
the importance of rotation in this regard (Meynet & Maeder 2000; Maeder & Meynet 2000; 
Heger & Langer 2000; Heger, Langer, & Woosley 2000). Yet these models, untested as they 
may be, provide the only means for linking a variety of observational studies to astrophysical 
interpretation, such as the determination of the initial mass function (Massey 1998) and the 
determination of turn-off masses in clusters containing massive stars (Massey, Waterhouse, 
& DeGioia-Eastwood 2000, Massey, DeGioia-Eastwood, & Waterhouse 2001). 

Herrcro ct al. (1992) first called attention to a significant mass discrepancy between 
the masses derived from modern stellar atmosphere models, and that inferred from stellar 
evolution codes, for the highest mass stars, in the sense that the evolutionary tracks predict 
a mass as much as 2 times larger. (See also Herrero, Puis, & Villamariz 2000.) Burkholder, 
Massey, & Morrell (1997) recently tried to resolve this mass discrepancy by using data on 
massive spectroscopic binaries. They found good agreement between the binary masses and 
the evolutionary tracks up to 25^1©. Some higher mass systems did show significant lower 
masses than the evolutionary tracks would suggest, but these systems were either at or near 
their Roche lobes, suggesting that significant mass might have been lost from the system. 
(See also Penny, Gies, & Bagnuolo 1999.) 

Massey & Hunter (1998) recently obtained spectra of 65 of the bluest, most luminous 
stars in the R136 cluster located at the heart of the 30 Doradus nebula in the LMC. The 
majority of these stars proved to be of spectral type 03, the hottest, most luminous, and 
most massive stars known. Four of the stars showed widely spaced, double absorption lines, 
indicative of spectroscopic binaries caught at favorable phases. Since the evolutionary masses 
of these stars were very high, it was thought that these would be excellent candidates for 
additional study, in the hopes that orbit solutions would resolve the mass discrepancy once 
and for all. 
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2. Observations and Reductions 
2.1. The Data 

The data were all obtained as part of a Cycle 8 HST program, GO-8217. We utilized 
30 orbits, organized into 15 "visits" of 2 consecutive orbits, with each visit separated by 
carefully planned intervals. We used the Space Telescope Imaging Spectrograph (STIS) in 
imaging mode for photometry of our binaries, and with a medium-dispersion grating for 
spectroscopy of each binary. 

Each visit began with a pair of short images (total integration time of 2.2 seconds) 
obtained with the "long-pass" filter (which cuts off all light < 5500A) and centered on the 
middle of the cluster. The 28 arcsec by 51 arcsec field of view always contained our four 
binaries, but changes in the spacecraft roll angle resulted in some of the out-lying R136 stars 
either being included or not. The spatial scale was 0.05 arcsecs per pixel. 

After centering up on a fairly isolated offset star (Melnick 34 = R136-8; see Fig. 1 
of Massey & Hunter 1998), the telescope was precisely offset to each of the four binaries 
in turn, in order to obtain spectra with the G430M grating centered at A4451. These 
exposures covered the wavelength rang e 4310A to 4590A at a dispersion of 0.28A per pixel, 
with a nominal 1.5-pixel resolution of 0.42A (28 km s~^). This wavelength region was 
selected to include as many spectral lines as possible in a single wavelength setting; i.e., 
H7, He I A4471, and He H A4542. We used the 0.2 arcsec by 52 arcsec slit. Wavelength 
cahbration was obtained for each exposure using the default iJST/STIS scheme, resulting in 
a new calibration exposure prior to the first spectroscopic exposure for each orbit. 

Near the end of each visit, another pair of short-exposure images (total integration time 
2.6 seconds) were obtained to continue the photometric monitoring. The total elapsed time 
for each visit was 3.2 hours. 

In order to maximize our phase coverage for all periods of interest, we adopted a clever 
scheme developed by Abi Saha for observing Cepheids as part of the distance-scale key 
project: we designed our program so that each visit was separated in time according to a 
geometrical progression, corresponding to 

gap(n)[days] = 0.5 * 1.175", n = 1, 14 

The multiplicative and geometrical factors are chosen to assure good phase coverage for 
the shortest and longest periods of interest, given a predetermined number of observations. 
Simulations showed that this would provide uniformly good phase coverage for periods from 
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< 1 day through 30 days. The longest interval between successive visits was 3.7 days, and 
the shortest interval was 0.5 days. The observations spanned an observing "season" of 24.4 
days. In arranging the order of our visits, we put all the shorter intervals near the middle 
of the 24 day sequence, with progressively longer intervals to either end. The result of this 
is that any eclipses would tend to be found in the middle of the sequence, as we will find in 
Sec. 3.6 is indeed the case. 



2.2. Photometry 

For the STIS images, we used the "standard pipeline" fiat-fielded images prior to cosmic- 
ray (CR) rejection. They had been observed in "CR-split" mode, resulting in a pair of 
images, each of 1.1-1.3 sec duration. We chose to do photometry of all 60 images, and 
rejected CR events by comparing the photometry from each of the CR-split images. If the 
photometry agreed to within the photometric errors (3(j), we averaged the results; if not, 
we assumed that the smaller magnitude was spurious due to a CR event. The photometry 
was measured using IRAF's aperture photometry routine, with a radius of 2 pixels and sky 
being determined from the modal value within an annulus lying 5-10 pixels from the star's 
center. The point-spread-function has a full- width at half-radius of 1.7 pixels, and we found 
that the centers were well determined by a Gaussian fit to the radial profile. 

In order to provide more accurate differential photometry for our binaries, and to search 
for other light variables, we performed photometry of 59 of the brightest, most isolated stars 
in the cluster, including our four binaries. Of the 59 stars, one happened to lie outside the 
ficld-of-view for two of the visits, due to the slightly changing roll-angle of the spacecraft 
over the course of the observations. After eliminating possible variables (based upon whether 
agreement frame-to- frame was within the photometric errors), we found that there were 
indeed photometric zeropoint changes of ~ 0.04 mag from frame-to-frame, with the drifts 
being secular rather than random in time: i.e., over the course of several days the zero- 
point would drift upwards and then downwards again. This is consistent with the claim 
that "instrumental [in] stability" limits absolute photometry with STIS to 5% (Table 16.3 
in Leitherer et al. 2001). In any event, our use of STIS as an N-star photometer (with 
N~ 58 — 59) resulted in photometry with high precision for temporal changes. 

We have applied a single adjustment to the photometric zero-point to have the photom- 
etry roughly match V; a single value worked well as the colors of these stars are all very 
similar (Hunter et al. 1997). 
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2.3. Spectroscopy 

Owing to the extreme crowding in the R136 cluster, the spectra for each of the four 
binaries had to be re-extracted with a small aperture centered on each target, with care 
being taken to select "clean" sky on either side of the slit. We actually reduced our spectra 
two ways. For one version, we began with the pipeline final two-dimensional spectrum 
(flat-fielded, CR-removed, geometrically-corrected, and wavelength cahbrated) and simply 
re-extracted the spectra of our binaries. For the other, we began with the fiat-fielded version, 
and extracted our spectra using the usual IRAF routines; for these we also applied our own 
dispersion solution determined from the wavelength-calibration frames. We found in practice 
that our own reductions provided better signal-to-noise. 

In measuring the spectra, we fit double lines using two Vogt functions; this provided 
excellent fits to even blended double lines, and avoids issues with pair-blendings, always a 
concern with orbit solutions of broad-lined stars (see discussion in Burkholder et al. 1997). 
The single-lined phases were measured with a single Vogt function for consistency. Both 
versions of the spectra were measured independently, and the results compared; no systematic 
differences were found, and the comparison proved useful mainly for eliminating any spurious 
measurements due to noise spikes. Our measuring error was 5-10 km s~^. 

3. Results 

We provide the photometric data and radial velocity data in Tables 1 and 2. In this 
section we describe what we learned from the spectra of each binary, and provide details of 
the orbit solutions and analysis of the light-curve information. We also report the discovery 
of five more light-variables in the R136 cluster. We begin by explaining our methods. 

3.1. Methodology 

For the period searches, we used a version of the Lafier & Kinman (1965) routine, which 
relies upon point-to-point smoothness after the data (either photometric or radial velocity) 
are sorted in phase according to a trial period. For the orbit solutions, we initially solved 
for each component independently using a modified version of the differential corrections 
program of Wolfe, Horak, & Storer (1967), fixing only the period. In all cases we found 
solutions consistent with circular orbits, which we expect given the short periods and high 
masses (see Sect. 4). We then determined the best values of the orbital semi-amplitudes K 
"center-of-mass" 7-velocities by running a non-linear least-square routine based upon the 
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grid-search program of Bevington's (1969), with the eccentricity, time of conjunction, and 
periods fixed for the two components. This then maximizes the precision of our determination 
of 7 and K. As a reminder, we expect that the 7-velocities of the two components of a massive 
binary will differ due to the photospheric outflow velocities of the stellar winds (Massey & 
Conti 1977). We characterize the agreement between the orbit and velocity data by Rl, 
computed from the goodness-of-fit x- 

The spectral types and magnitude difference between the components can be determined 
from the best double-lined phases. Our spectral types differ slightly from those of Massey & 
Hunter (1998), as these STIS spectra have considerably higher resolution and better signal- 
to-noise. The relative strengths of He II A4542 and He I A4471 helped establish the spectral 
class; the absolute visual magnitudes were consistent with all of these stars being dwarfs. 
The magnitude difference could be measured from the relative fluxes of the H7 hues at 
double-lined phases, as the equivalent width of H7 is fairly insensitive to T^s for dwarfs this 
hot. 

We measured the projected rotation velocity vsvai for each component by using the 
appropriate model atmosphere lines and convolving these with rotational velocities until we 
obtained the best fit to a fine. The model atmosphere code is that described by Kudritzki & 
Puis (2000), and we used models computed by PM for stars of similar spectral types in the 
LMC Given the intrinsic line widths, we found we could not measure rotational velocities 
smaller than 90-100 km s^^. In the following, we will compare the rotational velocities to 
that expected on the basis of synchronous rotation, computed using the stellar radii and 
orbital period. 

We thus know a great deal about the physical parameters of these stars, which is helpful 
in reducing the number of free parameters in interpreting the light-curves. We adopt a 
distance modulus (m — Af)„ = 18.5, in accord with Westerlund (1997) and van den Bergh 
(2000). Accurate values of the reddening were determined by Massey & Hunter (1998) 
based upon the mulit-band HST photometry of Hunter et al. (1997). Thus the absolute 
visual magnitudes are known (Table 1 of Massey & Hunter 1998) for the combined systems, 
and, combined with the magnitude differences measured from the spectra, tells us My for 
each component in the binary. The spectral type allows us to assign effective temperatures; 
we adopt the scale of Chlebowski & Garmany (1991). (We will discuss this in more detail 
in Sec. 4.) The effective temperature then determines the bolometric correction (Vacca, 
Garmany, & ShuU 1996), and hence we know the bolometric luminosity of each star. 

We used the fight curve synthesis code GENSYN (Mochnacki & Doughty 1972) to 
produce model light curves. Our approach was to make a constrained fit using as much data 
as possible from the spectroscopic results as described above. The orbital parameters were 
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taken from the spectroscopic solution, and the physical parameters were estimated from 
the spectral classifications of the stars. We set the stellar temperatures according to the 
spectral classification. Wc then estimated the physical fluxes and limb darkening coefficients 
from tables in Kurucz (1979) flux models and Claret (2000), respectively. The observed flux 
ratios together with the adopted effective temperatures yield estimates of the ratio of stellar 
radii. Then, for a given input value of the polar radius of the primary Rp, we calculated 
the secondary radius. Each trial run of GENSYN was set by two independent parameters, 
the system inclination i and primary polar radius. For each run, we attempt to match two 
observables: the absolute visual magnitude of the system and the eclipse depths and widths. 
The best fit solutions are those with the calculated My of the system that also matched the 
eclipse depths. Models that fit the eclipse depths could be made for a range of inclinations, 
however the My for such models would greatly diverge from that calculated using the well 
known distance and reddening to the LMC. The quoted errors in inclination derive from an 
estimated error of ± 0.15 mag on our My values, which is based upon the uncertainty in 
the LMC distance modulus (van den Bergh 2000) and a modest error in the photometry and 
reddening (Massey & Hunter 1998). 

For all four systems, the stars are well contained within their Roche surfaces. In the 
case of R136-39, no eclipses are seen. However, the crucial phases where we might expect 
eclipses lack observations. Therefore, we only quote an upper limit on the inclination. Any 
inclination above this value would result in eclipses that would be both too deep and wide 
to agree with the current observations. 

As a further check on the models, we independently estimated the orbital inclination 
simply from geometry, after measuring the eclipse depths. For this, we adopt a modest 
correction for limb darkening using a linear coefficient (Al-Naimiy 1978; van Hamme 1993), 
but ignore refiection and other second-order effects. Given the poor sampling in phase-space 
we expected only modest agreement with the models, but in fact the agreement was excellent, 
giving us high confidence that the orbital inclinations are well determined. 

We estimate the errors on the physical parameters, including a 10% uncertainty in the 
effective temperature scale for 0-type stars (Conti 1988). Nevertheless, our parameters are 
well determined, in large part because the inverse dependence of the stellar radii on effective 
temperature is partially canceled by the dependence of the bolometric correction on effective 
temperature, and in fact the uncertainties in Am dominate the errors on the stellar radii. 

L/Lq = {R/Ref X (Teff/To)^ 

R/R^ = (L/L^f' X {Tjnsf (1) 
For solar-type stars, a small error in estimating the effective temperature Tes would result 
in a large error in R/Rq, since there is such a steep inverse relationship with Teff- However, 
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for hot stars we must make a very large correction for the bolometric correction (BC) in 
determining L/ Lq starting with My, and that the BC has a very steep dependence on T^s 
as well: 

BC = 27.66 - 6.84 X log Teff, 

from Vacca et al. (1996). Thus substituting this into equation (1), and adopting Tq — 
5770°K, we find 

R/Rq = 871.67 X IQ-^y/^ X 
By propagation of errors we find then that 
4/i?o = <ff X (-550.90 X X T-fl-^'^^f + al,^ x (-401.42 x IQ-^v/^ x Tj-^^'^f 



3.2. R136-38 

Six of the 15 spectra showed double lines at H7 and He 11 A4542. There is a large mag- 
nitude difference between the two components, which we measure as Am=1.0±0.2 mag. The 
He IA4471 line clearly follows the motion of the secondary, but never showed any component 
due to the primary. Thus we have a good measurement of the motion of the secondary 
for both single-lined and double-lined phases. Occasionally the He I hne was too noisy to 
measure. We find that the single-lined phases follow the motion of the primary very well, 
but we will not give those any weight in the orbit solutions. The primary is of spectral type 
03 V, while the secondary is of spectral type 06 V, as judged by the relative strengths of 
He 1 and He 11 during the double-lined phases. The photometry shows well pronounced dips 
of ~ 0.2 mag indicative of eclipses. 

Period searches on the radial velocity data (both primary and secondary) and on the 
photometry all revealed the same period, 3.39 days. We present the orbital parameters in 
Table 3, and show the radial velocity curves and orbit fit in Fig 1. We find that the rotational 
velocities are consistent with synchronous rotation. Analysis of the light-curve finds a well- 
determined orbital inclination oii — 79° ± 1°. We show the light curve data, and the model 
fit, in Fig 1 as well. 

This is a very interesting system, containing stars of extremely high masses, with the 
03 V primary having a mass of blM.Q. This is higher than the mass determined from any 
other binary, exceeding the mass of even that of Plaskett's star (51A40, Bagnuolo, Gies, 
& Wiggs 1992).^ We compare these to the masses derived from the evolutionary tracks in 
Table 3, and find excellent agreement. 



^There are two other contenders for the recorder-holder. (1) HD 92740 is a WN7 Wolf-Rayet star in 
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3.3. 



R136-39 



This system consists of an 03 V + 05.5 V pair. Of the 15 spectra, 11 show double 
hnes at H7 and He II A4542, and it is clear the the He I A4471 comes purely from the 
secondary. The magnitude difference between the two components is 0.45 ±0.1 mag. The 
period is 4.06 days, as determined from the radial-velocities. Inspection of the unphased 
photometry did not show any obvious signs of an eclipse, but when phased according to the 
radial velocity information we find that there is poor coverage near the important phases. 
This allows us only to assign an upper limit on i, and our modeling suggests that the orbital 
inclination must be less than 75° to account for the lack of deeper eclipses. This places 
only lower limits on the masses on the system, of 27A^0 and 21A4.Q. The line widths are 
consistent with synchronous rotation. 



This system consists of an 03 V + 03 V pair. Of the 15 spectra, 11 showed double lines 
at H7 and He II A4542, and there was no trace of He I in any of our spectra. The magnitude 
difference between the two components is modest, Am=0.2 ± 0.1. The light-curve shows a 
very deep eclipse (0.5 mag), suggesting that we are viewing this system at a very favorable 
inclination. 

A period search of both the radial velocity and photometry data yielded the same value, 
P = 2.89 days. We give the orbital parameters in Table 5, and show the orbit and light- 
curve in Fig. 3. We find i = 85.4°. The masses are among the highest seen in a spectroscopic 
binary, 40^1© and 33 Ai®. The line widths are consistent with synchronous rotation. 



The last system we discuss consists of two 05.5 V + 05.5 V stars with equal brightness. 
Of the 15 spectra, 12 showed double lines at H7, He I A4471, and He II A4542. However, 



which the absorption hnes and emission hnes follow the same orbital motion. By combining spectra from six 
He I lines, Schweickhardt et al. (1999) report the detection of a secondary absorption spectrum, which they 
describe as 09 III. The double-lined orbit solution yields a mass of 55.3 ± I.SMq for the Wolf-Rayet star. 
(2) HD 93205 is an 03 V + 08 V pair recently studied by Antokhina et al. (2000). The mass of the 03 V 
primary was found to be 32-154A^0, with a most probable value of 45 Mq. We are indebted to Doug Gies 
for reminding us of these systems. 



3.4. 



R136-42 



3.5. 



R136-77 
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since we could not distinguish the primary from the secondary, it was necessary to rely upon 
the velocity differences between the two components in order to find the period; a period 
search of the photometry yielded identical results, and strong (0.4 mag) eclipses are evident. 
The period was then used to phase the radial velocity data, and assign velocities to one star 
or the other. We designate the slightly more massive star the primary. The radial velocity 
data, orbit solution, and light-curve is shown in Fig. 4. The line widths are again found to 
be consistent with synchronous rotation. 

3.6. Other Light Variables 

In doing the photometry, we found five additional stars with significant photometric 
variations. We show the light-curves in Fig. 5. R136-07 (Mclnick 39), R136-15 (Melnick 30), 
R136-24 (R136a7), and R136-25 all show signatures of what might well be eclipses. The star 
R136-08 (Melnick 34) shows a much more puzzling behavior, changing by several tenths of a 
magnitude over the course of 3 weeks. The variations appear to be periodic. The spectrum 
of R136-08 mimics that of a Wolf-Rayet star, although Massey & Hunter (1998) argue that it 
is simply a "super Of star whose very high luminosity and stellar winds result in a spectrum 
dominated by emission. Spectroscopic and photometric monitoring of all five of these stars 
has been proposed for Cycle 11 with HST. 

4. Discussion 

Our analysis of the four R136 binaries have revealed orbital masses that are among the 
highest ever directly measured via this simple application of Kepler's second law. We can 
use these new data to compare to the masses derived from the evolutionary tracks^. We have 
included these values in Tables 3-6, and present the H-R diagrams (HRDs) in Fig 6. The 
agreement between the masses derived from these tracks, and the actual measured masses, 
is excellent for the three cases with eclipses. 

In placing the stars in the H-R diagram, we have chosen to adopt the effective temper- 
ature scale of Chlebowski & Garmany (1991). This is similar to the scale given by Conti 
(1988), who provides a critical discussion, and concludes that the absolute (but not relative) 



''We have used the older Schaerer et al. (1993) evolutionary tracks as these were the last set made public 
by the Geneva group that includes normal (rather than enhanced) mass-loss rates. Newer models including 
rotation are becoming available, but as yet none with the metallicity appropriate to the LMC. However, we 
expect that this effect will be small near the ZAMS, as suggested by Fig. 6 in Maeder & Meynet (2001). 
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uncertainties in the scale are about 10%. Since that time, more modern atmospheric models 
have been used to analyze a number of 0-type stars in the Milky Way, LMC, and SMC; see, 
for example. Puis et al. (1996). Such studies led Vacca et al. (1996) to propose a new effective 
temperature scale, which is ~6% hotter than the Conti (1988) calibration, and ~3% hotter 
than the Chlebowski & Garmany (1991) scale for the spectral types discussed in the current 
paper. We note that there are no spectral type to effective temperature scales determined for 
stars in the LMC and SMC, and that Vacca et al. (1996) restricted their study to Galactic 
stars. Thus refinements in the effective temperature scale will change the location of the 
stars in the HRD but probably the error bars in Fig. 6 are realistic. Note that we have only 
included the uncertainty in My in estimating the errors on the evolutionary track masses in 
Tables 3-6; were we to include the uncertainty in the effective temperature scale as well, the 
percentage error would roughly double. 

The stars do fall shghtly to the left of where we expect in the HRDs. Massey & Hunter 
(1998) found ages of 1-2 Myr for the R136 cluster, with the larger value corresponding 
to the cooler effective temperature scale, which we have used here. Yet the components 
in all four of our binaries lie on or near the ZAMS, to higher effective temperatures than 
the 2 Myr isochrones shown in Fig. 6. We do not have a ready explanation for this. We 
note that all four of these systems are relatively close pairs. Comparison of the orbital 
separations asini with the stellar radii (both quantities appear in Tables 3-6) reveals that 
the components are typically separated by ~ 2x the sum of the radii. This is sufficiently 
close that tidal forces have played a significant role. We see ready evidence of this in that 
there must have been some dynamical evolution for the orbits to be circular and the stars 
to be locked in synchronous rotation.^ Such tidal forces may have affected the evolution of 
the stars (providing additional heating on the envelopes), and in that case, these systems 
are not telling us as much as we would like about single stars. It is hoped that some of the 
stars identified in Sect. 3.6 may provide example of massive binaries with longer periods. 

Such systems would also help determine if the new generation of rotating stellar models 
do better than the standard non-rotating models. Macdcr & Meynet (2000) have recently 
invoked rotation to explain the discrepancy in masses between the evolutionary tracks and 
stellar atmosphere calculations; they note that tracks which fail to include rotation may 
overestimate the mass by as much as 50% in cases of high rotation and luminosity. The 
short periods of our binaries have resulted in slow rotation due to tidal forces, leading to 
little difference in the masses predicted by non-rotating and rotating models. 



^The time expected for tidal forces to circularize an orbit can be estimated using equation (2) of Shu & 
Lubow (1981); adopting parameters appropriate to high-mass stars with convective cores (Zahn 1975, 1977) 
leads to ~ 500, 000 years for these systems. 
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Fig. 1. — Radial velocity and light curves for R136-38. In the radial velocity curve, the 
solid points designate the primary; the asterisks denote the secondary, and the open circles 
denote the measurements at single-lined phases, which were not included in the orbit fit. In 
the plot of the light-curve, we include the the la error bars in the photometry, although for 
R136-38 these are comparable to the size of the points. The solid lines shows the results of 
our best model fit. 




Fig. 2. — Radial velocity and light curves for R136-39; symbols are the same as in Fig. 1. 
The model light-curve fit is for the maximum possible eclipse depth consistent with the data; 
we use this for setting an upper limit on the orbital inclination. 




Fig. 3.- 



Radial velocity and light curves for Rl 36-42; symbols are the same as in Fig. 1. 
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Fig. 4.- 



Radial velocity and light curves for Rl 36-77; symbols are the same as in Fig. 1. 
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Fig. 5. — Light-curves for five otlier suspected binaries. R136-08 (Melnick 34) show gradual 
variations with a period of 20 days; the others show changes typical of eclipses. 
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Fig. 6. — The locations of each star in the H-R diagram is shown by by the large filled circles. 
Error bars denote the change in location due to a 5% error in the effective temperature scale, 
and due to uncertainties in magnitude difference in the two components. The evolutionary 
tracks (sohd fines) are for z — 0.008 and come from Schaerer et al. (1993). The dotted fines 
are isochrones for 1, 2, and 3 Myr, computed from the same models. 
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Table 1. Journal of Observations, Radial Velocities 



Image Modified Orbital Radial Velocities (km s ^) 
Name Julian Date Phase'^ Primary Secondary Single 



R136-38 



O59001020 


51622.80 


0.41 




548 


202 


O59002020 


51626.29 


0.44 




473 


219 


O59003020 


51629.63 


0.43 




486 


219 


O59004020 


51632.17 


0.18 


118 


626 




O59005020 


51633.40 


0.54 






300 


O59006020 


51634.03 


0.73 


450 


-118^ 




O59007020 


51634.55 


0.88 




-32 


403 


O59008020 


51635.39 


0.13 


145 


581 




O59009020 


51635.97 


0.30 


112 


650 




O59010020 


51636.91 


0.58 






308 


O59011020 


51638.31 


0.99 




242 


279 


O59012020 


51639.32 


0.29 


117 


647 




O59013020 


51641.25 


0.86 


418 


-96 




O59014020 


51643.45 


0.51 






265 


O59015020 


51647.15 


0.60 




19 


332 


R136-39 


O59001030 


51622.84 


0.56 




190 


262 


O59002030 


51626.34 


0.42 


14P 


412 




O59003030 


51629.67 


0.24 


80 


514 




O59004030 


51632.20 


0.86 


429 


53 




O59005030 


51633.42 


0.16 


102 


481 




O59006030 


51634.04 


0.31 


97 


509 




O59007030 


51634.56 


0.44 




343 


218 


O59008030 


51635.42 


0.65 


429 


57 




O59009030 


51635.98 


0.79 


454 


9 
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Table 1 — Continued 



Image 


Modified 


Orbital 


Radial Velocities (km s ^) 


Name 


Julian Date 


Phase'' Primary 


Secondary 


Single 


O59010030 


51636.92 


0.02 




288^ 


251 


O59011030 


51638.36 


0.38 


120 


446 




O59012030 


51639.36 


0.62 


431 


69 




O59013030 


51641.31 


0.10 


148 


448 




O59014030 


51643.46 


0.63 


425 


45 




O59015030 


51647.21 


0.56 


355^ 


148^ 




R136-42 


O59001040 


51622.86 


0.27 


2 


578 




O59002040 


51626.34 


0.48 






250 


O59003040 


51629.69 


0.64 


468 


15 




O59004040 


51632.20 


0.51 






279 


O59005040 


51633.45 


0.94 






290 


O59006040 


51634.05 


0.15 


57 


548 




O59007040 


51634.53 


0.31 


18 


586 




O59008040 


51635.43 


0.62 


470 


-7 




O59009040 


51635.99 


0.82 


529 


-29 




O59010040 


51636.93 


0.14 


56 


535 




O59011040 


51638.38 


0.64 


505 


-2 




O59012040 


51639.38 


0.99 






274 


O59013040 


51641.39 


0.69 


530 


-32 




O59014040 


51643.48 


0.41 


128 


496 




O59015040 


51647.22 


0.70 


552 


-73 




R136-77 


O59001050 


51622.92 


0.29 


-33 


615 




O59002050 


51626.41 


0.15 


87 


561 




O59003050 


51629.74 


0.92 


454 


77 
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Table 1 — Continued 



Image 


Modified 


Orbital 


Radial Velocities (km s ^) 


Name 


Julian Date 


Phase'' 


Primary 


Secondary 


Single 


O59004050 


51632.23 


0.25 


-30 


602 




O59005050 


51633.47 


0.90 


458 


116 




O59006050 


51634.07 


0.22 


-11 


600 




O59007050 


51634.63 


0.52 






271 


O59008050 


51635.49 


0.98 






275 


O59009050 


51636.03 


0.27 


-36 


623 




O59010050 


51636.97 


0.77 


560 


-78 




O59011050 


51638.44 


0.55 






276 


O59012050 


51639.43 


0.08 


99 


467 




O59013050 


51641.37 


0.11 


49 


516 




O59014050 


51643.53 


0.26 


-45 


593 




O59015050 


51647.27 


0.25 


-56 


607 





^Given half weight in orbit solution. 

''Based upon the orbital parameters given subsequently. 



Table 2. Journal of Observations, Photometry 



Image Name 


MJD 


r(,ioD-C)o 


rllC)D-oy 


R136-42'= 


R136-77^ 


O59001010 


51622, 


.78 


14.39 


14.61 


14, 


.61 


15 


.27 


O59001060 


51622, 


.93 


14.36 


14.61 


14, 


.60 


15 


.30 


O59002010 


51626, 


.28 


14.35 


14.62 


14, 


.83 


15 


.40 


O59002060 


51626, 


.47 


14.51 


14.60 


14, 


.94 


15 


.28 


O59003010 


51629, 


.61 


14.34 


14.60 


14, 


.61 


15 


.28 


O59003060 


51629, 


.76 


14.39 


14.59 


14, 


.62 


15 


.29 


O59004010 


51636, 


.89 


14.38 


14.63 


14, 


.61 


15 


.24 


O59004060 


51636, 


.98 


14.34 


14.61 


14, 


.60 


15 


.22 


O59005010 


51638, 


.30 


14.54 


14.61 


14, 


.60 


15 


.60 


O59005060 


51638, 


.46 


14.47 


14.60 


14, 


.60 


15 


.42 


O59006010 


51639, 


.30 


14.32 


14.59 


14, 


.87 


15 


.70 


O59006060 


51639, 


.46 


14.32 


14.60 


15, 


.01 


15 


.25 


O59007010 


51641, 


.24 


14.35 


14.59 


14, 


.60 


15 


.47 


O59007060 


51641, 


.40 


14.34 


14.62 


14, 


.61 


15 


.27 


O59008010 


51643, 


.41 


14.47 


14.61 


14, 


.62 


15 


.24 


O59008060 


51643, 


.55 


14.49 


14.60 


14, 


.65 


15, 


.26 


O59009010 


51647, 


.14 


14.36 


14.60 


14, 


.59 


15 


.22 


O59009060 


51647, 


.30 


14.32 


14.60 


14, 


.61 


15 


.26 


O59010010 


51632, 


.16 


14.35 


14.61 


15, 


.12 


15 


.29 


O59010060 


51632, 


.26 


14.33 


14.60 


14, 


.87 


15 


.26 


O59011010 


51633, 


.39 


14.49 


14.61 


14, 


.60 


15 


.27 


O59011060 


51633 


.51 


14.36 


14.60 


14, 


.86 


15 


.33 


O59012010 


51634, 


.02 


14.34 


14.60 


14, 


.61 


15 


.24 


O59012060 


51634, 


.10 


14.32 


14.62 


14, 


.59 


15 


.24 


O59013010 


51634, 


.52 


14.33 


14.60 


14, 


.60 


15 


.49 


O59013060 


51634, 


.64 


14.34 


14.61 


14, 


.63 


15 


.60 


O59014010 


51635, 


.36 


14.34 


14.61 


14, 


.62 


15 


.28 


O59014060 


51635, 


.51 


14.33 


14.61 


14, 


.59 


15 


.68 


O59015010 


51635, 


.96 


14.35 


14.59 


14, 


.58 


15 


.25 


O59015060 


51636, 


.04 


14.33 


14.61 


14, 


.61 


15 


.23 



'^The photometric error for R366-38 is la — 0.007 mag. 
•^The photometric error for R366-39 is la — 0.009 mag. 

■^The photometric error for R366-42 is la = 0.010 mag. 
■^The photometric error for R366-77 is la = 0.027 mag. 
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Table 3. Orbit Solutions and Physical Parameters R136-38 



Parameter 



System 



Primary 



Secondary 



Orbital 
P (days) 



p. conjunction 



(MJD) 



7 (km s ^) 

K (km s-^) 

mp/ms 

Rl (km s-i) 

asini (i?o) 

msin^i {Mq) 

Spectral and Photometric 

Spectral Type 

Tefr(°K) 

Am 

Mv 

BC 

Radius (/?„) 

^sync 

vsini 

Echpse depths 

inclination (geometry) [deg] 

inclination (GENSYN model) [deg] 

Masses 
m (Mq) orbit 

m {M.q) evolutionary tracks 



3.39 (adopted) 
0.00 (adopted) 
51621.40 ±0.03 



0.41 



40.2 



1.0 ±0.2 
-5.3 



79.0 ± 1.0 
79.0 ±1.0 



278.2 ± 0.4 
174.7 ±0.5 

2.3 
11.7 
53.8 ±0.2 



272.3 ± 0.3 
424.8 ±0.2 

20.3 
28.5 
22.1 ±0.1 



03 V 06 V 

48, 500 ± 4, 850^ 42, 200 ± 4, 220^ 



-4.9 ±0.2"^ 
-4.4 ± 0.3^ 
-9.3 ± 0.4^^'^ 
9.3 ± LO^'*" 
110 ± 12 
130 ± 20 
0.23 



-3.9 ±0.2"^ 
-4.0 ±0.3^^ 
-7.9 ± 0.4^''' 
6.4 ±0.7 
76 ±9 
90 ±20 
0.20 



56.9 ±0.6 
53 ±5" 



23.4 ± 0.2 
29 ±2^ 



'^Adopting a 10% uncertainity in the spectral-type to Teg relationship. 
Errors (anti-) correlated. 
^Errors on the masses from the evolutionary tracks are based solely on the errors in My- 
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Table 4. Orbit Solutions and Physical Parameters R136-39 



Parameter 



System 



Primary 



Secondary 



Orbital 
P (days) 



p. conjunction 



(MJD) 



7 (km s ^) 

K (km s-^) 

mp/ms 

Rl (km s-i) 

asini (i?o) 

msin^i {Mq) 

Spectral and Photometric 

Spectral Type 

Tefr(°K) 

Am 

Mv 

BC 

Radius (/?„) 

^sync 

T;sini 

Echpse depths 
inclination (geometry) [deg] 
inclination (GENSYN model) [deg] 

Masses 
m (Mq) orbit 

m {M.q) evolutionary tracks 



4.06 (adopted) 
0.00 (adopted) 
51620.59 ±0.01 



0.754 ± 0.002 



37.5 



0.45 ±0.1 
-5.2 



< 72 



0? 



< 75.0 



271.8 ±0.3 
200.8 ±0.5 

7.1 
16.1 
24.5 ±0.1 



262.0 ± 0.3 
266.3 ±0.4 

9.5 
21.4 
18.5 ±0.1 



03 V 05.5 V 

48, 500 ± 4, 850^ 43, 200 ± 4, 320^ 



-4.7 ±0.1'^ 
-4.4 ± 0.3^ 
-9.0 ± 0.3^'^ 
8.1 ± 0.6^''' 
83 ±6 
< 100 
< 0.05? 



-4.2 ±0.1'^ 
-4.1 ±0.3^ 
-8.3 ± 0.3^^'^ 
7.1 ±0.6 
71 ±6 
< 100 



> 27.2 

46 ±2^ 



> 20.5 
34 ±2^ 



'^Adopting a 10% uncertainity in the spectral-type to Teg relationship. 
Errors (anti-) correlated. 
^Errors on the masses from the evolutionary tracks are based solely on the errors in My- 
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Table 5. Orbit Solutions and Physical Parameters R136-42 



Parameter 



System 



Primary 



Secondary 



Orbital 
P (days) 



p. conjunction 



(MJD) 



7 (km s ^) 

K (km s-^) 

mp/ms 

Rl (km s-i) 

asini (i?o) 

msin^i {Mq) 

Spectral and Photometric 

Spectral Type 

Tefr(°K) 

Am 

Mv 

BC 

Radius (/?„) 

^sync 

vsini 

Echpse depths 

inclination (geometry) [deg] 

inclination (GENSYN model) [deg] 

Masses 
m (Mq) orbit 

m {M.q) evolutionary tracks 



2.89 (adopted) 
0.00 (adopted) 
51622.07 ±0.01 



0.809 ±0.001 



35.5 



0.2 ±0.1 
-5.1 



85.5 ± 1.0 
85.4 ±0.5 



276.0 ± 0.3 
278.2 ± 0.4 

4.7 
15.9 
39.9 ±0.1 



268.0 ± 0.3 
343.7 ±0.4 

12.8 
19.6 
32.3 ±0.1 



03 V 03 V 

48, 500 ± 4, 850^ 48, 500 ± 4, 850^ 



-4.4 ±0.1"^ 
-4.4 ± 0.3^ 
-8.8 ± 0.3^'^ 
7.4 ± 0.8^'^ 
102 ±9 
100 ± 20 
> 0.40 



-4.2 ±0.1"^ 
-4.4 ± 0.3^ 
-8.6 ± 0.3^^'^ 
6.7 ±0.7 
93 ±7 
130 ± 30 
0.52 



40.3 ±0.1 
42 ±2^ 



32.6 ±0.1 
39 ±3" 



'^Adopting a 10% uncertainity in the spectral-type to Teg relationship. 
Errors (anti-) correlated. 
^Errors on the masses from the evolutionary tracks are based solely on the errors in My- 
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Table 6. Orbit Solutions and Physical Parameters R136-77 



Parameter 



System 



Primary 



Secondary 



Orbital 
P (days) 



p. conjunction 



(MJD) 



7 (km s ^) 

K (km s-^) 

mp/ms 

Rl (km s-i) 

asini (i?o) 

msin^i {Mq) 

Spectral and Photometric 

Spectral Type 

Tefr(°K) 

Am 

Mv 

BC 

Radius (/?„) 

^sync 

vsini 

Echpse depths 

inclination (geometry) [deg] 

inclination (GENSYN model) [deg] 

Masses 
m (Mq) orbit 

m {M.q) evolutionary tracks 



1.88 (adopted) 
0.00 (adopted) 
51622.37 ±0.01 



0.904 ± 0.001 



24.1 



0.0 ±0.2 
-4.5 



83.0 ± 1.0 
83.0 ±1.0 



273.8 ± 0.3 
309.2 ± 0.4 

17.5 
11.5 
28.3 ±0.1 



276.6 ± 0.3 
341.9 ±0.3 

14.8 
12.7 
25.6 ±0.1 



05.5 V 05.5 V 

43, 200 ± 4, 320^ 43, 200 ± 4, 320^ 



-3.7 ±0.2"^ 
-4.1 ±0.3^ 
-7.8 ± 0.4^^'^ 
5.8 ± 0.5^'^ 
124 ± 13 
140 ± 20 
0.45 



-3.7 ±0.2"^ 
-4.1 ±0.3^ 
-7.8 ± 0.4^'^ 
5.8 ±0.5 
124 ± 13 
130 ± 20 
> 0.35 



28.9 ±0.3 
28 ± r 



26.2 ±0.3 
28±0.r 



'^Adopting a 10% uncertainity in the spectral-type to Teg relationship. 
Errors (anti-) correlated. 
^Errors on the masses from the evolutionary tracks are based solely on the errors in My- 



